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Thermodynamics relies on the possibility to describe systems composed of a large number of 
constituents in terms of few macroscopic variables. Its foundations are rooted into the paradigm of 
statistical mechanics, where thermal properties originate from averaging procedures which smoothen 
out local details. While undoubtedly successful, elegant and formally correct, this approach carries 
over an operational problem: what is the precision at which such variables are inferred, when 
technical/practical limitations restrict our capabilities to local probing? Here we introduce the 
local quantum thermal susceptibility, a quantifier for the best achievable accuracy for temperature 
estimation via local measurements. Our method relies on basic concepts of quantum estimation 
theory, providing an operative strategy to address the local thermal response of arbitrary quantum 
systems at equilibrium. At low temperatures it highlights the local distinguishability of the ground 
state from the excited sub-manifolds, thus providing a method to locate quantum phase transitions. 


The measurement of temperature is a key aspect in 
science, technology, and in our daily life. Many inge¬ 
nious solutions have been designed to approach different 
situations and required accuracies [T]. What is the ul¬ 
timate limit to the precision at which the temperature 
of a macroscopic state can be determined? An elegant 
answer to this question is offered by estimation theory n- 
13]: The precision is related to the heat capacity of the 
system [51 [Sj. 

In view of the groundbreaking potentialities offered 
by present-day nanotechnologies iiHni and the need to 
control the temperature at the nano-scale, it is highly 
relevant to question whether the heat capacity is still 
the relevant (fundamental) precision limit to small-scale 
thermometry. The extensivity of the heat capacity is a 
consequence of the growing volume-to-surface ratio with 
the size m- However, at a microscopic level such con¬ 
struction may present some limitations [m |T5| . More¬ 
over a series of theoretical efforts recently concentrated 
on a self-consistent generalization of the classical ther¬ 
modynamics to small-scale physics, where quantum ef¬ 
fects become predominant [TBHUj. In particular, a lot of 
attention has been devoted to the search for novel meth¬ 
ods of precision nanothermometry that could exploit the 
essence of quantum correlations [251l28j . In this context, 
the possibility to correctly define the thermodynamical 
limit, and therefore the existence of the temperature in 
the quantum regime, has been thoroughly investigated. 
It has been shown that the minimal subset of an in¬ 
teracting quantum system which can be described as a 
canonical ensemble, with the same temperature of the 
global system, depends not only on the strength of the 
correlations within the system, but also on the tempera¬ 
ture itself |29H31j . Using a quantum information-oriented 
point of view, this phenomenon has been also highlighted 
in Gaussian fermionic and bosonic states, by exploiting 
quantum fidelity as the figure of merit [32| [33]. Fur¬ 
thermore, the signihcant role played by quantum corre¬ 
lations has been recently discussed with specihc attention 


to spin- and fermonic-lattice systems with short-range in¬ 
teractions [34] . 

In this paper we propose a quantum-metrology ap¬ 
proach to thermometry, through the analysis of the lo¬ 
cal sensitivity of generic quantum systems to their global 
temperature. Our strategy does not assume any con¬ 
straint neither on the structure of the local quantum 
state, nor on the presence of strong quantum fluctuations 
within the system itself, while it rather moves from the 
observation that the temperature is a parameter which 
can be addressed only via indirect measurements. Specif¬ 
ically we introduce a new quantity which we dub Lo¬ 
cal Quantum Thermal Susceptibility (LQTS), and which 
gauges how efficiently the thermal equilibrium of a sys¬ 
tem is perceived by its subparts. More precisely, given 
a quantum system AB in a thermal equilibrium state, 
the LQTS ©_4 is a response functional which quantifies 
the highest achievable accuracy for estimating the system 
temperature T through local measurements performed on 
a selected subsystem A of AB (see Fig. [^. 

The LQTS is in general not extensive with respect to 



FIG. 1: The operationally grounded strategy associ¬ 
ated to the LQTS functional. A composite quantum sys¬ 
tem S is in thermal equilibrium with a bath at temperature 
T. The LQTS functional 0 measures the highest achievable 
accuracy in the estimation T under the hypothesis to perform 
only local measurements on the subsystem A of AB. 
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the size of A, yet it is an increasing function of the latter, 
and it reduces to the system heat capacity in the limit 
where the probed part coincides with the whole system 
AB. In the low-temperature limit, we shall also see that 
the LQTS is sensitive to the local distinguishability be¬ 
tween the ground state and the first excited subspace 
of the composite system Hamiltonian. In this regime, 
even for a tiny size of the probed subsystem, our func¬ 
tional is able to predict the behaviour of the heat ca¬ 
pacity and in particular to reveal the presence of criti¬ 
cal regions. This naturally suggests the interpretation of 
as a sort of mesoscopic version of the heat capacity 
which replaces the latter in those regimes where exten- 
sivity breaks down. 


RESULTS 


The functional. Let us consider a bipartite quan¬ 
tum system AB at thermal equilibrium, composed of two 
subsystems A and B, and described by the canonical 
Gibbs ensemble pp = /Zp. Here H = Hj\ + Hq + 
is the system Hamiltonian, which in the general 
case will include both local (ie., iL _4 and H^) and in¬ 
teraction (i.e., terms, while Zp — Tr[e“^^] = 

Si denotes the associated partition function (/? = 

l/kgT is the inverse temperature of the system, kg the 
Boltzmann constant, and {Ei} the eigenvalues of H). In 
this scenario we are interested in characterizing how the 
actual temperature T is perceived locally on A. 

For this purpose we resort to quantum metrology [35] 
and define the LQTS of subsystem A as 


&a[pp\ '■= 8 lini 


1 - E{pj, 
^2 


( 1 ) 


where F (p, a) = yTp ayTp] is the fidelity between 

two generic quantum states p and a [36l [37] . The quan¬ 
tity 0 corresponds to the quantum Fisher information 
(QFI) [31 m for the estimation of /3, computed on the re¬ 
duced state p'p = Trglpp]. It gauges how modifications 
on the global system temperature are affecting A, the 
larger being &j\[pp\ the more sensitive being the sub¬ 
system response. More precisely, through the quantum 
Cramer-Rao inequality, &a[Pi 3 \ quantifies the ultimate 
precision limit to estimate the temperature T, by means 
of any possible local (quantum) measurement on sub¬ 
system A. In the specific, it defines an asymptotically 
achievable lower bound. 


AT-^>A:BTVyiV6^M, (2) 

on the root-mean-square error AT'^ = y^E[(T®®* — T)^] 
of a generic local estimation strategy, where T®®* is the 
estimated value of T, E[a;] is the expectation value for 
a random variable x, and N is the number of times the 
local measurement is repeated. 


By construction, &a[pp\ is a positive quantity which 
diminishes as the size of A is reduced, the smaller being 
the portion of the system we have access to, the worse 
being the accuracy we can achieve. More precisely, given 
A' a proper subset of A, we have &A'[pp\ A &a[pp\- In 
particular, when A coincides with the whole system AB^ 
equation Q reaches its maximum value and becomes 
equal to the variance of the energy, 

<Sab[pp\ = Ev[ppH^] - Tv[ppHf , (3) 


which depends only on the spectral properties of the sys¬ 
tem and which coincides with the system heat capac¬ 
ity [3113] (note that, rigorously speaking, the LQTS quan¬ 
tifies the sensitivity of the system to its inverse tempera¬ 
ture /3; the corresponding susceptibility to T = I/(A:b/ 3) 
gets a kg^T~‘^ correction term, which also enters the 
standard definition of the heat capacity). 

An explicit evaluation of the limit in equation 0 can 
be obtained via the Uhlmann’s theorem [35] (see Methods 
for details). A convenient way to express the final result 
can be obtained by introducing an ancillary system A'B' 
isomorphic to AB and the purification oi pp defined as 

e-/3SA2 

- -7^^\Ei)AB ® \Ei)A'B' , (4) 

where H = J2iEi\Ei)AB{Ei \ is the spectral decomposition 
of the system Hamiltonian. It can then be proved that 

&A[pp] = &AB[pp] - E > (5) 

j<k ^ 


\pp) = E 


where {|e^)} are the eigenvectors of the reduced density 
matrix Tr^/[jp^) (p^|] living on ABB', obtained by tak¬ 
ing the partial trace of \pp) with respect to the ancillary 
system A', while {Xj} are the corresponding eigenvalues 
(which, by construction, coincide with the eigenvalues 
of pf). 

Equation 0 makes it explicit the ordering between 
&a[pp\ and &ab[pp\'- the latter is always greater than 
the former due to the negativity of the second contri¬ 
bution appearing on the rhs. Furthermore, if H does 
not include interaction terms {i.e., = 0), one can 

easily verify that &a[pp] reduces to the variance of the 
local Hamiltonian of A, and is given by the heat capac¬ 
ity of the Gibbs state jZ^ which, in this spe¬ 
cial case, represents p-p, i.e., &a[pp\ = — 

Tr[p^FGl]^ = --^InZ-p. Finally we observe that in 
the high-temperature regime (/3 —?> 0) the expression 0 
simplifies yielding 


&a[pp\ = jjAt: 


H\ - 2PHa 


TrB[g^ 

d-B 




■ o{p^ 


( 6 ) 

where cIa and djg denote the Hilbert space dimensions of 
A and B respectively, and we defined Ha = 
having set, without loss of generality, Tr[77] = 0. 
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FIG. 2: Schematic representation of the low-energy 
spectrum for a generic many-body quantum system. 

For simplicity the ground-state (gs) energy Eq is set to zero. 
Here Hi denotes the normalized projector on the eigenspaces 
of energy Ei, which can be rii-fold degenerate. 


A measure of state distinguishability. In the 
low-temperature regime, the LQTS can be used to char¬ 
acterize how much the ground state of the system AB 
differs from the first excited subspaces when observing 
it locally on A. This is a direct consequence of the fact 
that the QFI (which we used to define our functional) 
accounts for the degree of statistical distinguishability 
between two quantum states (in our case the reduced 
density matrices and differing by an infinites¬ 

imal change in the investigated parameter (in our case 
the inverse temperature /3). Therefore for /3 —> oo, the 
LQTS can be thought as a quantifier of the local distin¬ 
guishability among the lowest energy levels in which the 
system is frozen. 

To clarify this point, let us consider the general sce¬ 
nario depicted in Fig.[^ where we only discuss the physics 
of the ground state (with energy Eq = 0 ) and of the low¬ 
est excited levels with energy Ei bounded by twice the 
energy of the first excited level, Ei < 2Ei. The degener¬ 
acy of each considered energy eigenstate is denoted by rii. 
From equation ([^ it then follows that up to first order 
in the parameter we get 

&a[P0] = E (l-Tr[Po^n^]) +0 . 

no V / 

I 

(7) 

Here H^ = Trgpi], where H^ is the normalized pro¬ 
jector on the degenerate subspace of energy Ei. More¬ 
over, is the span of the local subspace associated 

A. 

to the ground state, ie., 

A 

Hg' = '^^^iPj\4>j)A{4‘j\: ^-nd njf being the number of 
non-zero eigenvalues {pj > 0 ) of Hjf. 

Equation Q can be interpreted as follows. Our ca¬ 
pability of measuring /3 relies on the distinguishability 
between the states p^ and Pp+g, with e (3. In the zero- 
temperature limit, the system lies in the ground state and 
locally reads as Hj)^, while at small temperatures the low¬ 
est energy levels start to get populated. If their reduced 
projectors Hd (j > I) are not completely contained in the 


span of Hq^, that is Tr[PQ^n^] 7 ^ 1 , there exist some lo¬ 
cal states whose population is null for T = 0 and greater 
than zero at infinitesimal temperatures. Such difference 
implies that the first order in ©_4 does not vanish. On 
the contrary, if the reduced projectors H^ are completely 
contained in the span of Hq^, that is Tr[Po^Uf] = 1 , we 
can distinguish from Pp+g only thanks to infinitesi¬ 
mal corrections (!I(exp(—2/3i?i)) to the finite-valued pop¬ 
ulations of the lowest energy levels (see Methods for an 
explicit evaluation of the latter). In conclusion, the quan¬ 
tity 6 yi[pp->.oo] acts as a thermodynamical indicator of 
the degree of distinguishability between the ground-state 
eigenspace and the lowest energy levels in the system 
Hamiltonian. 

LQTS and phase estimation. A rather stimulat¬ 
ing way to interpret equation ([^ comes from the observa¬ 
tion that, in the extended scenario where we have purified 
AB as in equation Q , the global variance ([^ formally co¬ 
incides with the QFI SabA'B'Hp^^^)) associated with the 
estimation of a phase p which, for given /3, has been im¬ 
printed into the system ABA'B' by a unitary transforma¬ 
tion with PP being the analogous of El on the 

ancillary system A'B', i.e., 6 ^b[pp] = ^ABA'B'{\pf'^)) 

where |pp^^) = e“®'^'‘^/^|pp) [H [351 I3S]- Interestingly 
enough a similar connection can be also established with 
the second term appearing in the rhs of equation ( 1 ^: 
indeed the latter coincides with the QFI 5 'b^'B'(| pp^Q) 
which defines the Cramer-Rao bound for the estimation 
of the phase ip of |pp^^), under the constraint of having ac¬ 
cess only on the subsystem BA'B' (i.e., that part of the 
global system which is complementary to A). Accord¬ 
ingly we can thus express the LQTS as the difference be¬ 
tween these two QFI phase estimation terms, the global 
one vs. the local one or, by a simple rearrangement of the 
various contributions, construct the following identity 

&a[pp] + SbA'B'Hp^^^)) = , (8) 

which establishes a complementarity relation between the 
temperature estimation on A and the phase estimation 
on its complementary counterpart BA'B', by forcing their 
corresponding accuracies to sum up to the energy vari¬ 
ance (AiL^)p of the global system (|^. 

Local thermometry in many-body systems. We 
have tested the behaviour of our functional on two models 
of quantum spin chains, with a low-energy physics char¬ 
acterized by the emergence of quantum phase transitions 
(QPTs) belonging to various universality classes [30] . 

Specifically, we consider the quantum spin-1/2 Ising 
and Heisenberg chains, in a transverse magnetic field h 
and with a z-axis anisotropy A respectively, 

+ ( 9 ) 

i 

i7xxz = j 5 ][(«+i + afaf+,) + Aa/a/+,].( 10 ) 
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FIG. 3: The LQTS in the Ising and the Heisenberg models at low temperature. We numerically computed the LQTS 
of equation in the low-temperature limit for a chain with L = 12 sites in the following two cases: (a) the Ising model as a 
function of the transverse field h; (b) the Heisenberg XXZ chain as a function of the anisotropy A. The uppermost (red) curve 
corresponds to the global quantum thermal susceptibility, that is the heat capacity. The other curves stand for different sizes 
riA of the measured subsystem A of AB {ua increases along the direction of the arrow). The inset in panel (b) magnifies the 
data around A = — 1. In the XXZ model, the LQTS with ha ~ 1 can be proved to rigorously vanish. The inverse temperature 
has been fixed in both cases at /3 = 9. 


Here ct“ denotes the usual Pauli matrices {a = x, y, z) 
on the ith site, and periodic boundary conditions have 
been assumed. We set J = 1 as the system’s energy 
scale. At zero temperature, the model presents a Z 2 - 
symmetry breaking phase transition at | /ic | = 1 belonging 


to the Ising universality class. The Hamiltonian (10 1 has 


a critical behaviour for —1 < A < 1, while it presents 
a ferromagnetic or antiferromagnetic ordering elsewhere. 
In the latter case, the system exhibits a first-order QPT 
in correspondence to the ferromagnetic point Af = —1, 
and a continuous QPT of the Kosterlitz-Thouless type at 
the antiferromagnetic point Aaf = 1 . 

Figure displays the small-temperature limit of 
&a[pp\ for the two models above, numerically computed 
by exploiting expression (21) in Methods. We first ob¬ 


serve that, as expected, for all the values of h and A the 
LQTS monotonically increases with increasing the num¬ 
ber UA of contiguous spins belonging to the measured 
subsystem A. More interestingly, we find that even when 
A reduces to two or three sites, its thermal behaviour 
qualitatively reproduces the same features of the global 
system (represented in both models by the uppermost 
curve). In particular, even at finite temperatures and for 
systems composed of twelve sites, the LQTS is sensitive 
to the presence of critical regions where quantum fluc¬ 
tuations overwhelm thermal ones. The reminiscence of 


QPTs at finite temperatures has been already discussed 
via a quantum-metrology approach, through the analy¬ 
sis of the Bures metric tensor in the parameter space 
associated to the temperature and the external parame¬ 
ters m- The diagonal element of such tensor referring 
to infinitesimal variations in temperature, corresponds to 
the thermal susceptibility of the whole system. The latter 
quantity has been recently studied for the XY model [28] , 
showing its sensitivity to critical points of Ising univer¬ 
sality class. 

In the low-temperature regime such global sensitivity 
can be understood within the Landau-Zener (LZ) for¬ 
malism |32|- This consists of a two-level system, whose 
energy gap AE varies with respect to an external control 
parameter T, and presents a minimum AAmin in cor¬ 
respondence to some specific value Tc. Conversely, the 
global heat capacity ([^ may exhibit a maximum or a 
local minimum at Tc, according to whether AAmi„ is 
greater or lower than the value of AE* at which the 
expression &ab[p/3] = p is maximum in AA, re¬ 

spectively. Indeed it can be shown that {AEl"^)p for a 
two-level system exhibits a non-monotonic behaviour as 
a function of AE, at fixed (3 (see Methods). Quite re¬ 
cently, an analogous mechanisms has been also pointed 
out for the global heat capacity in the Lipkin-Meshkov- 
Glick model m- The LZ formalism represents a simpli- 
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fied picture of the mechanism underlying QPTs in many- 
body systems. However, by definition the temperature 
triggers the level statistics and the equilibrium properties 
of physical systems. Therefore, both the heat capacity of 
the global system and the LQTS of its subsystems 
are expected to be extremely sensitive to the presence of 
critical regions in the Hamiltonian parameter space. 

We also performed a finite-size scaling analysis of 
&a[pp\ as a function of the size of the measured sub¬ 
system. For slightly interacting systems, one expects the 
LQTS to be well approximated by the heat capacity of A 
(at least when this subsystem is large enough). The lat¬ 
ter quantity should scale linearly with its size n_ 4 . This 
is indeed the case for the Ising model ([^, where a direct 
calculation of (AiL^)^ can be easily performed [28]. Our 
data for the scaling of the stationary points of &ji\pp\ 
close to QPTs suggest that significant deviations from a 
linear growth with are present (see the Supplementary 
Information). This indicates that correlations cannot be 
neglected for the sizes and the systems considered here. 


DISCUSSION 

We have proposed a theoretical approach to tempera¬ 
ture locality based on quantum estimation theory. Our 
method deals with the construction of the local quantum 
thermal susceptibility, which operationally highlights the 
degree at which the thermal equilibrium of the global 
system is perceived locally, avoiding any additional hy¬ 
pothesis on the local structure of the system. This func¬ 
tional corresponds to the highest achievable accuracy up 
to which it is possible to recover the system temperature 
at thermal equilibrium via local measurements. Let us 
remark that, even if in principle the Cramer-Rao bound 
is achievable, from a practical perspective it represents 
a quite demanding scenario, as it requires the precise 
knowledge of the Hamiltonian, the possibility to iden¬ 
tify and perform the optimal measurements on its sub¬ 
systems, and eventually a large number of copies of the 
system. However, in this manuscript we have adopted a 
more theoretical perspective, and focused on the geomet¬ 
rical structure of the quantum statistical model underly¬ 
ing local thermalization. 

In the low-temperature regime our functional admits 
an interpretation as a measure of the local state distin- 
guishability between the spaces spanned by the Hamil¬ 
tonian ground state and its first energy levels. Further¬ 
more we established a complementarity relation between 
the highest achievable accuracy in the local estimation of 
temperature and of a global phase, by showing that the 
corresponding accuracies associated with complementar¬ 
ity subsystems sum up to heat capacity of the global 
system. Finally, we considered two prototypical many- 
body systems featuring quantum phase transitions, and 
studied their thermal response at low temperatures. On 


one hand, we found that optimal measurements on local 
systems provide reliable predictions on the global heat 
capacity. On the other hand, our functional is sensitive 
to the presence of critical regions, even though the total 
system may reduce to a dozen of components and the 
measured subsystem to one or two sites. Let us remark 
that most of the results presented herewith do not refer to 
any specific choice of the interaction Hamiltonian, 
between A and B. As an interesting implementation of 
our scheme, we foresee the case of non-thermalizing in¬ 
teractions HSHm, whose potentialities for precision ther¬ 
mometry have been already unveiled. 

We conclude by noticing that, while in this article we 
focused on temperature, the presented approach can be 
extended to other thermodynamic quantities (like en¬ 
tropy, pressure, chemical potential etc.). Furthermore 
it seems plausible to adopt quantum estimation strate¬ 
gies to tackle the problem of providing self-consistent 
definition of heat and work for microscopic systems (see 
Ref. [IS] and references therein). One of the main dif¬ 
ficulties to this end derives indeed from the fact that 
these thermodynamic functionals are processes rather 
than state variables: accordingly there are no quan¬ 
tum observables associated with them. In this respect 
a quantum-estimation strategy not explicitly referred to 
a specific quantum observable, but rather bearing the 
geometrical traits of the Hilbert space associated to the 
explored systems, appears as a valid chance to close this 
gap. 


METHODS 

Derivation of equation Let us recall the def¬ 

inition of the LQTS for a given subsystem A of a global 
system S at thermal equilibrium: 

6 ^M= 81 imi—( 11 ) 

£—>■0 s 

where T (p, tr) = yTp ayTp] is the fidelity between 

two generic quantum states p and a. According to the 
Uhlmann’s theorem [38], we can compute J- as 

T {pj,p-^+,) = max |(p^|p^+J|, (12) 

|p^>.|P^+e> 

where the maximization involves all the possible purifica¬ 
tions \p-p) and Ipp+g) of p-p and Pp+g respectively through 
an ancillary system a. A convenient choice is to set 
a = BA'B', with A'B' isomorphic to AB. We then 
observe that, by construction, the vector |pp) of equa¬ 
tion besides being a purification of pp, is also a par¬ 
ticular purihcation of p^. We can now express the most 
generic purihcation of the latter as 

P-/3W2 

(13) 


I/Op) = (lyt o F)|pp) = (1.4 o V) 
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where V belongs to the set of unitary transformations 
on a, where represents the identity operator on the 
system X, and where in the last equality we introduced 
the vector jlij) = Y.i\Ei) ab\E^)a'B', \Ei)AB being the 
eigenvectors of H. We can thus write the fidelity (12) as 


X = 


l/3+e 


max 

V 




1^'B')] \Pl^) 


(14) 

Since we are interested in the small-e limit, without loss 
of generality we set V = exp(ien), with fl being an 
Hermitian operator on the ancillary system a. It comes 
out that, up to corrections of order 0{e^), the LQTS 
reads 


&a[pp] = &Ab[pp] +4iiun|Tr[p^n^] - Tr[p^n]^ 

+ l^{pp\P^H]\pp)^ , (15) 


p'p. To the latter scenario belongs the case in which 
H = Hj^+Hi 3 , where one can easily prove that the LQTS 
reduces to the variance of the local Hamiltonian LQi, i.e. 
6 ^ = Ti[ppH^] — Tr[p^iL_ 4 ]^ (notice that the non-zero 
eigenvalues of p^ are Xi = e~P^' jZ-^ which correspond 
to |e,) = \Et) ® |pf), being Hj, = 

= Tr[e“^'^^] and |p®) the purification of p® through 
the ancillary system B'). The expression above can be 
also rewritten as 

GaIpp] = GabIpp] - E ’ (20) 

j<k ^ 

which can be cast into equation ([^ by simply exploiting 
the fact that the system is symmetric with respect to the 
exchange of AB with A'B'. 

It is finally useful to observe that LQTS can be also 
expressed in terms of the eigenvectors of A, = 
as: 


where we have defined pp = Tr_ 4 [|pp)(pp|]. By differ¬ 
entiating the trace with respect to XI we determine the 
minimization condition for it, yielding 

(fl-a;)p^-bp^(H-w) = Q, (16) 


= ( 21 ) 

J,k 

where we have used the Schmidt decomposition of |pp), 
with respect to bipartition Aa, 


with u! = Tr[p^fl] and Q = —^[H',Pp], H' being the 
analogous of H which acts on A'B' (by construction 
H\pp) = H'\pp)). Equation explicitly implies that 
XI does not depend on w, which, without loss of general¬ 
ity, can be set to zero. Moreover, it enables to rewrite 
the LQTS in equation lH) as 

Ga[pp] = Gab[pp] - 4Tr[p^ H"] . (17) 


The solution of the operatorial equation (16) can be 
found by applying the Lemma presented at the end of 
this section, yielding 


n = flo + ^ {PH'R + RH'P) 


i 



' + 00 

^(-l)"(Pp^®"iL'p^" 


_n—l 



(18) 


\Pp) — ^ ( \/^\9i)\^i) ■ 


( 22 ) 


In particular, expression (21) can be exploited in order to 


numerically compute the LQTS, for instance when deal¬ 
ing with quantum many-body systems (see Fig. and 
the discussion in the Supplementary Information). 


Lemma. For any assigned operators X, Y satisfying the 
eguation 

XW + WX = V , (23) 

admit solutions of the form 

W = Wo + X^^YR + RYX'^^ 

+ 00 

+ ^(-i)"X’"py(x®("+b - R ), (24) 

n—0 


with Oq being an operator which anti-commutes with XI, 
poGm being the Moore-Penrose pseudoinverse of pp to 
the power m, while R being the projector on kernel of p^, 
and P = Ha — R being is complementary counterpart. By 
substituting this expression in equation we finally 
get 


where is the Moore Penrose pseudoinverse of X, 

R is the projector on kernel of X, P = t — R (^1 indi¬ 
cates the identity matrix) and Wq is an generic operator 
which anti-commutes with Y. Furthermore if X and Y 
are Fdermitian, equation (231 admits solutions which are 
Hermitian too: the latter can be expressed as 


&a[pp\ = 2E - ^APpHf , (19) 

i.fc ^ 

where p^ = i® ^be spectral decomposition 

of p‘p, sharing the same spectrum with p^. The expres¬ 
sion above holds for both invertible and not invertible 


W = Wn-hX^^YR + RYX^^ 


(25) 


1 


+ 00 


^(-i)"x”py(x®("+i) -R) + h.c. 


^ n—0 


where now Wq is an arbitrary Hermitian operator which 
anti-commutes with Y. 
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Proof. Since (23) is a linear equation, a generic solution 
can be expressed as the sum of a particular solution plus 
a solution Wq of the associated homogeneous equation, 
i.e. an operator which anti-commute with X, 


X= 0 . 


(26) 


A particular solution W of equation (231 can be always 
decomposed as 

W = RWR + RWP + PWR + PWP. (27) 

Notice that by definition RX = XR = O, where O iden¬ 
tifies the null operator. Multiplying (231 on both sides by 
R, one gets the condition RYR = O. The operator W, 
solution of equation (23), is defined up to its projection 
on the kernel subspace, that is 

W' = W + RWR XW' + W'X = Y . (28) 

Therefore, without loss of generality we can set 

RWR = O. (29) 


Multiplying equation (23) by on the right side and 
repeating the same operation on the left side we get: 


RWP = RYX^\ 
PWR = X®^YR. 


(30) 

(31) 


On the other hand, PWP satisfies the relation 

PWP = PYX^^ - PYR - X{PWP)X^^. (32) 

This equation can be solved recursively in PWP and 
gives 


PWP = ^(-l)"A”Pr(A®(i+") - R), 


(33) 


n=0 


thus concluding the first part of the proof. The second 
part of the proof follows simply by observing that, if X 
and Y are Hermitian, and if IT solves equation (23), then 
also its adjoint counterpart does. Therefore for each so¬ 
lution W of the problem we can construct an Hermitian 
one by simply taking (IT -|- W^)/2. □ 

Second-order term corrections to LQTS. In 

the low-temperature regime (/3 —)■ oo), we have computed 
the second-order correction term to the LQTS, that is of 
6.4 [P/ 3 ] in equation ([^, and found: 


'Einf\ _ 
no 


-2/3£i 


-2-bTr[Po^nf-b nf 

+00 


nj^(i + Po^)] 


-2^(-l)"Tr 


n—0 






with Ek > El and where the series in n is meant to con¬ 


verge to 1/2 when Tr 


n+1 


= 1 


, z.e., 



FIG. 4: The two-level Landau-Zener model. A sketch of 
the behaviour of the two eigenenergies (Ti, E 2 ) as a function 
of some control parameter F. The gap AE = E 2 —E 1 displays 
a pronounced minimum in correspondence of a given F'^ value. 


a;” = 5 - In order to van¬ 
ish, this second-order correction term requires a stronger 
condition with respect to one necessary to nullify the 
first-order term in the LQTS, equation ([^. It is given 
by = Hg^, and corresponds to the requirement that 
the system ground state must be locally indistinguishable 
from the first excited level. 

Heat capacity in the 2-level Landau-Zener 
scheme. Here we discuss the simplified case in which 
only the ground state (with energy Eq) and the first ex¬ 
cited level (with energy Ei) of the global system Hamil¬ 
tonian H play a role. In particular, we are interested 
in addressing a situation where the ground-state energy 
gap AE = El — Eq may become very small, as a function 
of some external control parameter T {e.g., the magnetic 
field or the system anisotropy). A sketch is depicted in 
Fig. El. and refers to the so called Landau-Zener (LZ) 
model [H]. This resembles the usual scenario when a 
given many-body system is adiabatically driven, at zero 
temperature, across a quantum phase transition point. 

In correspondence of some critical value Tc, the gap 
is minimum. For a typical quantum many-body system, 
such minimum value AAmin tends to close at the ther¬ 
modynamic limit and a quantum phase transition occurs 
(notice that F^ may depend on the system size). Here¬ 
after, without loss of generality, we will assume Eq = 0 
and take Ei — AE so that the system heat capacity ^ 
reduces to: 


^ ^ ^ ni{AE)^e-f^^^ (niAEe-^^^\ 

- no + nie-/5AiJ " Vno+nie-Z^AS ) ' 

(34) 

Here no and ni are the degeneracy indexes associated to 
the levels Eq and Ei, respectively. Notice that &ab[pp\ 
is always non-negative and exhibits a non-monotonic be¬ 
haviour as a function of AE, at fixed /3. Indeed it is im¬ 
mediate to see that &ab\pp\ —>■ 0 in both limits AE —>■ 0 
and AE —>■ -boo. For fixed /3, hq and ni, the heat capac¬ 
ity displays a maximum in correspondence of the solution 





















FIG. 5: Heat capacity in the LZ model. The emergence 
of extremal points in the behaviour of &ab[pp\ as a function 
of increasing control parameter F is associated to the relative 
size of AiJmin with respect to /S.E*. The red arrows denote the 
changing of Ai? during a typical LZ protocol. One realizes 
that, if Ai?min > Ai?* one peak will appear (left), while if 
AjSniin < AS* one peaks will appear (right). 

of the transcendental equation 

d&Aisipp] „ r, . , ^AB _ ^1 (2 +/3 AS) 

9AS "no(/3AS-2)- 

In particular, for ng = rii = 1 the latter relation is ful¬ 
filled for AS* « 2.3994//3, while for no = 2,ni = 1 it is 
fulfilled for AS* « 2.2278/(3. 

It turns out that the behaviour of the heat capacity 
as a function of increasing F in a two-level LZ scheme 
depends on the relative sizes of AS* and as 

pictorially shown in Fig. a) if ASmin > AS*, then 
QabIpp] will exhibit a maximum in correspondence of 
Fci b) if ASmin < AS*, a maximum at F/ corresponding 
to AS = AS* will appear, followed by a local minimum 
at Fc and eventually by another maximum at F 2 where 
the former condition occurs again. Since AS* is a func¬ 
tion of (3, and ASmin depends on the system size, the 
point of minimum gap can be signaled by a maximum 
or by a local minimum depending on the way the two 
limits L —>■ -boo (thermodynamic limit) and (3 —> -boo 
(zero-temperature limit) are performed. 
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SUPPLEMENTARY INFORMATION 

Scaling of the LQTS with the dimension 
of subsystem A. 


a power-law behaviour of the peak value at least for very 
small subsystems, scaling as 

&A[pp] - {n/Lr, (36) 


The local quantum thermal susceptibility (LQTS) for 
the subsystem ^ of a given system, as defined in equa¬ 
tion Q, is a quantity which is tricky to be evaluated 
numerically. Apart from the exponential growth of the 
Hilbert space, extrapolating the limit £ —0 typically re¬ 
quires high accuracies in the diagonalization procedure. 
This would generally limit the study of local thermome¬ 
try in the many-body context up to very small systems. 
The expression that we derived in equation (211 circum¬ 
vents the latter problem and enables an easier manipula¬ 
tion of &a[P/3]- However one still needs the full spectrum 
of the reduced density matrix since all its eigenvec¬ 
tors have to be contracted with the product of the global 
equilibrium state times the system Hamiltonian, ppH. 
This makes the whole analysis not straightforward, even 
for free-fermion systems as is the case for the Ising model 
^^ising- For this reason we resort to an exact diagonaliza¬ 
tion technique. 

Let us concentrate on the scaling of the LQTS with 
the size n _4 of the measured subsystem for the Ising spin 
chain The maximum of ©^[pp] as a function of 
and for different system sizes L is shown in Fig. As 
noticed in the main text within the LZ discussion, the 
peak is related to the appearance of the critical point, 
but deviations from its position at the thermodynamic 
limit (he = 1) are expected, due to the mutual interplay 
of finite-size and finite-temperature effects. To balance 
them, we choose a size-dependent temperature /3 = 3L/4. 
While with exact diagonalization we cannot go beyond 
the study of systems with size L = 12, our data suggest 



with an exponent a « 2. 

We performed the same analysis also for the XXZ 
chain (101, and focused on the behaviour of &a[pp] 
around the two critical points of the model, in corre¬ 
spondence to the two local minima. As we did above, 
we fix a size-dependent temperature /3 = 3L/4 and study 
systems of size up to L = 12. In the left panel of Fig. 
we concentrate on the ferromagnetic point at Af = —1. 
Similarly to the critical point of the Ising model, a power- 
law scaling of the type in equation (361 seems to emerge 
at small n^, with an exponent a « 3. The scaling anal¬ 
ysis at the antiferromagnetic point around Aaf = 1 is 
less clear and probably requires larger system sizes (right 
panel). We notice the appearance of a cusp-like feature 
at ua = L/2, which reflects the different behaviour of 
&a[pp\ for A « Aaf, depending on whether ua A L/2 or 
UA > L/2, as is clearly visible in the right panel of Fig. 
We have to stress that the departure of the two bunches 
of curves becomes more evident at large values of /3, while 
it tends to disappear when increasing the temperature. 



FIG. 7: Same quantity and parameters as in Fig. but for 
the XXZ chain. The data in the two panels are for the minima 
of the LQTS close to the critical points Af = — 1 (left) and 
Aaf = +1 (right). The dashed line in the left panel denotes 
the behaviour 6^ ~ (n/L)® and is plotted as a guideline. 


Analysis of the heat capacity by considering 
only the first excited levels. 


FIG. 6: Peak value of [p/j] for the Ising model, as a function 
of the dimension of the measured subsystem ua- Different 
data sets are for various system lengths L. The dashed line 
denotes a power-law behaviour Ga ~ (’I'/L)^, and is plotted 
as a guideline. The temperature is chosen as a function of the 
system size, as /3 = 3Z//4. 


The heat capacity of the global system is quantified by 
the variance of the energy ([^, and it can be calculated 
easier than the LQTS, since it does not involve any par¬ 
tial tracing over a portion of the system and depends only 
on the spectral properties. In particular, at low temper¬ 
atures (3 —>■ -boo, the relevant contributions to &ab[pp\ 
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FIG. 8: Heat capacity for the Ising chain, without breaking 
the Z 2 symmetry. The upper panel shows the hrst energy 
gaps AEj = Ej — i?o as a function of the transverse field h. 
The dashed green line denotes the gap AEi in the thermody¬ 
namic limit. The lower panel shows the global quantum ther¬ 
mal susceptibility &ab[pp\ for the full spectrum (red), and 
when taking into account the first eigenenergies only. Here 
we consider a chain with L = 12 and a temperature /3 = 9. 


will be provided only by the first low-lying energy levels, 
as discussed within the Landau-Zener context. Here we 
provide an analysis of the role of such excited states in 
the two many-body systems that we address in this work. 

We start from the Ising chain The low-lying spec¬ 
trum is shown in the upper panel of Fig. where we 
plot the gaps AEj — Ej — Eq between the ground-state 
energy Eq and those of the first excited states Ej, j > 0. 
Since we are considering the full Hilbert space of the sys¬ 
tem, for L —> 00 the ground state in the ferromagnetic 
side {\h\ < he) is doubly degenerate, while in the para¬ 
magnetic side {\h\ > he) a gap opens up monotonically 
as AEi = 2\h — he\- At finite values of L (data with 
symbols), the double degeneracy survives only at h = 0 
and Ai monotonically increases with h > 0. 

The scenario emerging by only keeping contributions 
to &ab[pp\ coming from the ground state and the first 
excited state, while neglecting any other excited level, is 
quite clear. In the thermodynamic limit and for finite 
/3, the heat capacity is rigorously zero for h < he and is 
finite for h > he, exhibiting a non-monotonic behaviour. 
The position h* of the maximum depends on the value of 
P and shifts toward he = 1 as long as /3 —>■ 00. Note how¬ 
ever that QabIpp] 0 . In summary, the position of 

the maximum depends on a competition between the fol¬ 
lowing two effects: i) decreasing L tends to shift the peak 
toward the ferromagnetic phase (h* < he); ii) decreasing 
P tends to shift the peak toward the paramagnetic phase 
{hi > he). 

In the bottom panel of Fig. we observe that, keep¬ 


ing only the ground and the hrst excited level in the 
computation of the energy variance (|^, already gives an 
excellent approximation to the exact value of &ab[pi3\. 
Therefore the above situation effectively applies at large 
P values. In the presence of a symmetry breaking mech¬ 
anism, one would hnd a situation analogous to the LZ 
scheme: the presence of a maximum in proximity of he 
would not be guaranteed, and a local minimum may ap¬ 
pear (according to the values of L and P). 

Let us now consider the XXZ-Heisenberg chain ( [To| . 
Analogously as above, the upper panel of Fig. [^displays 
the hrst energy levels in the full Hilbert space (we do 
not break the symmetry associated to the conservation of 
the global magnetization along the z axis). For A < Af, 
the ground state is fully polarized along z, and presents 
a double degeneracy at any length. At A = Aaf we 
see a cusp in the ground-state energy gap that is due 
to a level crossing (this feature persists at larger sizes 
L). This scenario characterizes the contribution to the 
heat capacity coming only by the ground state and the 
hrst excited state (black line in the bottom panel). In 
particular we see that, in this approximation, &ab[pp\ 
is zero for A < Af, while it becomes hnite at A > Af. 
On the other side, the Aaf is signaled by a cusp, which 
displays a minimum, and which raises as a consequence 
of the cusp in the ground-state energy gap. 

Contrary to the Ising model, if also other excited lev¬ 
els are taken into account, non-negligible corrections to 
the energy variance appear (see the bottom panel of 
Fig.§. In particular we note the emergence of a mini¬ 
mum around Af, due to the fact that &ab [pp] becomes h- 
nite also for < Af. The corner point at Aaf is smeared 
into a local minimum with a continuous derivative. The 
positions of the two minima are inhuenced by hnite-size 
and hnite-temperature effects. 
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FIG. 9: Same quantities and parameters as in Fig. but for 
the XXZ chain in the full many-body Hilbert space of L spins. 































